/* Integration points (quadratures)
   Jan Stebel
*/

#ifndef _QUADPOINT_H__
#define _QUADPOINT_H__


#include "intpoint.h"


/** Quadratic interpolation */
template<unsigned int dim>
class QuadPoint : public IntPoint<dim>
{
protected:
	using IntPoint<dim>::imesh;
	bool isLinear() { return false; };

public:
	QuadPoint() {};
	QuadPoint(Mesh *mesh, const Quadrature<dim> *q) { IntPoint<dim>::init(mesh, q, (dim+1)*(1+dim/2)); };
	void setGSF();
	void setSF();
	int mapDof(int elem, int dof, int iParent);
};









#endif
